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Abstract 



| We describe numerical calculation results for the probability distribution of the value of the 

monopole creation operator in the 577(2) lattice gluodynamics. We work in the maximal abelian 
projection. It occurs that at low temperatures, below the deconfinement phase transition, the 
■ maximum of the distribution is shifted from zero, which means that the effective constraint potential 

is of the Higgs type. Above the phase transition the minimum of the potential (the maximum of 
the monopole field distribution) is at the zero value of the monopole field. This fact confirms the 
existence of the abelian monopole condensate in the confinement phase of lattice gluodynamics, 
and agrees with the dual superconductor model of the confining vacuum. 

1 Introduction 

8* 

The monopole mechanism of the color confinement |2j, ^] is generally accepted by the 
lattice community. Still there are many open questions In the lattice gluodynam- 
ics it is very important to find the order parameter, constructed from the monopole field 
for the deconfinement phase transition. The first candidate is the value of the monopole 
condensate, which should be nonzero in the confinement phase and vanish at the phase 
transition. To study the monopole condensate, we need an explicit expression for the oper- 
ator Q m on(x), which creates the abelian monopole at the point x. The operator & m on(x), 
found for the compact electrodynamics with the Villain form of the action by Frohlich 
and Marchetti ||, was studied numerically in Q. In Section 2 we construct the analo- 
gous monopole creation operator for an arbitrary abelian projection of the lattice SU(2) 
gluodynamics. The numerical results presented in Section 3 are obtained for the maximal 
abelian projection. As shown by many numerical simulations for this projection the glu- 
odynamic vacuum behaves as the dual superconductor (see reviews || and references 
therein). In Q another form of the monopole creation operator was studied; it is shown 
that its expectation value vanishes in the deconfinement phase. The monopole creation 
operator, constructed in Section 2 is positive, and therefore its expectation value cannot 
vanish at T = T c . Still the results of the numerical studies of the effective potential for 
&mon presented in Section 3, clearly indicate that the monopole condensate exists in the 
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confinement phase of lattice gluodynamics and does not exist in the deconfinement phase. 
The analogous claim is done in ref.|R]], where the monopole condensate is calculated on 
the basis of the percolation properties of the monopole currents. The monopole creation 
operator in the monopole current representation is studied in ref.|ll|. Again, the monopole 



creation operator depends on the temperature as the disorder parameter. 

In Appendices A and B we prove that the operator used in the numerical calculations 
create the monopole in the lattice U(l) theory with the general form of the action. Ap- 
pendix C contains the brief description of the differential form notations on the lattice, 
these notations are used in Appendices A and B. 

2 Monopole Creation Operator 

First we give a formal construction of the monopole creation operator for the abelian 
projection of the SU(2) gluodynamics. Let us parametrize the SU (2) link matrix in the 
standard way: U% = cos^e*^; U% = sin^e^; U% = U%*; U% = -U™>; 
< < vr/2, — 7r < 6,x < 7T- The plaquette action in terms of the angles 0, 8 and x can 
be written as follows: 

S P = ^Tr U 1 U 2 U+U+ = S a + S n + S i , (1) 

where 

S a = COS0P COS 01 COS 02 COS 03 COS 04, (2) 

S n and S l describe the interaction of the fields 6 and x an d selfinteraction of the field 
X H^l; Op = 6 1 + 02 — 9s — 64, here the subscripts 1, ...,4 correspond to the links of the 
plaquette: 1 — > {x, x + /2}, 4 — > {x, x + z>}. For a fixed abelian projection, each term S a , 
S n and S l is invariant under the residual U(l) gauge transformations: 

Qxfi Qxfj, + ol x — a x+ f l , (3) 

Xxv Xxtx + Q-X + a x+ (x . (4) 

We define the operator which creates the monopoles at the point x of the dual lattice 
as follows: 

<$> mon (x) = exp {P[-S(9 P , ...) + S(6 P + W P (x), ...)]} , (5) 
the function Wp(x) being defined by eq. (|A.6[) . Substituting (|l|-|2|) into (|5|), we get 

<S>mon{x) =exp|^/3[-cos(0p) + cos(0p + VF P (x))]| , (6) 

where (5 = cos 0i cos 02 cos 03 cos 04 (3. Effectively the monopole creation operator shifts 
all abelian plaquette angles Op. 
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For the compact electrodynamics with the Villain type of the action the above definition 
coincides with the definition of Frohlich and Marchetti ||. For the general type of the 
action in compact electrodynamics we can use the described above construction. The 
proof is outlined in Appendices A,B. The gluodynamics in the abelian projection contains 
the compact gauge field 6 and the charged vector field x- The action (|l]) in terms of the 
fields 9 and x 1S rather nontrivial, and at the moment we can not prove that the above 
construction of the monopole creation operator is valid in this case. However, there is a 
proof for a similar Abelian - Higgs model, with the general type of the action, and this 
proof is analogous to one given in Appendix B. Moreover the numerical results, presented 
in the next section, clearly show that the introduced operator is the order parameter for 
the deconfinement phase transition. 



3 Numerical Results 

The numerical results are obtained on the lattice 4 • L 3 , for L = 8, 10, 12, 14, 16. The 



preliminary numerical results are published in ref. [13]. We extrapolate our results to the 
infinite volume since near the phase transition the finite volume effects are very strong. 
We impose anti-periodic boundary conditions in space directions for the abelian fields, 
since the construction of the operator f mon can be done only in the time slice with the 
anti-periodic boundary conditions. Periodic boundary conditions are forbidden due to the 



Gauss law: we input a magnetic charge into the finite box. Formally, equation (A. 5) for *D 
admits no solution in the finite box with periodic boundary conditions. To impose anti- 
periodic boundary conditions on the abelian fields the C-periodic boundary conditions 
should be imposed on the nonabelian gauge fields [14]. In the case of SU(2) gauge group 
the C-periodic boundary conditions are almost trivial: on the boundary we have U x a — > 
Q + U Xl u,ft>, Q = i<?2- To get the order parameter for the deconfinement phase transition we 
study the probability distribution of the operator & m0 n] we calculate the expectation value 
< 5($ — $>mon(x)) >■ The effective constraint potential, 

V eff m = -ln(< - Ij> mon(x )) >) (7) 

X 

has more physical meaning than the probability distribution. The calculation of V e ff(&) 
is time consuming, and we present our results for V(<J?), defined as follows: 



= -ln{< 6{<f> - $ m on{x)) >). (8) 

In Figs. 1(a) and 1(b) V(3>) is shown for the confinement and the deconfinement phases, 
the calculations being performed on the lattice 4 • 12 3 . 

In the confinement phase the minimum of V(<&) is shifted from zero, while in the 
deconfinement phase the minimum is at the zero value of the monopole field <£. We 
have used the positive operator, & m on(x) > @, however in the dual representation the 
creation operator of the monopole QA.5p is not positively definite: the sign is lost when we 
perform the inverse duality transformation (cf. eq. ( |A.5| ) with eq. ( |A.6j )). The potential 
shown in Fig. 1(a) corresponds to the Higgs type potential. The value of the monopole 
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field, $ c , at which the potential has a minimum is equal to the value of the monopole 
condensate. The potential shown in Fig. 1(b) corresponds to the trivial potential with a 
minimum at the zero value of the field: $ c = 0. The dependence of the minimum of the 
potential, <£ c , on the spatial size of the lattice, L, is shown in Fig. 2. The gauge fields 
are generated by the standard heat bath method. At each value of (3 2000 update sweeps 
are performed to thermalize the system. The maximal abelian projection Q corresponds 
to the maximization of the quantity R = Yl x ,ix Tr(JJ xlx <J3U£ a^). For our relatively small 
lattices the overrelaxation algorithm [|17], [l8| and simple local maximization method Q 
give approximately the same results, and we use the local maximization method. We stop 
our gauge fixing sweeps when Z = 10 -5 , here Z 18] is the lattice analogue of the quantity 
< \{d^ + iAfyA^l 2 > which should be zero if the maximal abelian projection is fixed 
exactly. We checked that the more precise gauge fixing do not change the results for <l? c 
inside the statistical errors. For each value of (3 at the lattice of definite size we use 100 
gauge field configurations separated by 300 Monte Carlo sweeps. For each configuration 
we calculated the value of the monopole creation operator at 20 randomly chosen lattice 
points. Therefore for each value of f3 at the lattice of definite size we have 2000 values of 
&mon- We use these values to calculate the quantity $ c (the maximum of the probability 
distribution p($ m on) is at <J> mori = $ c ). 

We fitted the data for <£ c by the formula <I> C = AL a + $>™f , where A, a and $* n ^ are 
the fitting parameters. It occurs that a = — 1 within statistical errors. Fig. 3 shows the 
dependence on (3 of the value of the monopole condensate, extrapolated to the infinite 
spatial volume, It is clearly seen that ffr 1 ™* vanishes at the point of the phase 

transition and it plays the role of the order parameter. 
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Appendix A 

Here we construct the monopole creation operator for the compact electrodynamics with 
the general type of the action. A similar construction exists for the compact Abelian-Higgs 
model with the general type of the action. First we perform the duality transformation of 
the partition function for the 4D lattice compact electrodynamics with the general type of 
the action S(dd + 2vr) = S(d<9), 
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V9exp{-S(d9)}, (A.l) 



We use the notations of the calculus of differential forms on the lattice |T9[ (see also 
Appendix C). The symbol J T>9 denotes the integral over all link variables 9. First consider 
the Fourier series for the Boltzmann factor: 



+ 7T +7T 

Z = const.Jv9 F(n)e^ n ' de \ F{n) = fvXexp{-S[X] -i(n,X)} . (A.2) 

— 7r n(c 2 )e^ 

— n 

Integrating over 9 we get the partition function of the dual theory: 

Z d = exp{-*S(d*k)}, (A.3) 

where n = 5k, *S(n) = — In F[n). We can represent Z as the following limit of the 
partition function for the Abelian-Higgs theory: 



+ 7T +00 

Z d = lkn J V*ip J V*B ^ exp{-*S{d* B /2tt) -r]\\*B -d*ip + 27r*k\\ 2 }, (A.4) 

here *S(d*B/2it) is the kinetic energy of the dual gauge field *B (the analogue of F^ u ) 
and the Higgs field exp{i*<^} carries magnetic charge, since it interacts via the covariant 



derivative with the dual gauge field *B. The Dirac operator [15], 



$>mon d {x) = e 1 ** ■ exp {-i(*D x , *B)} , S*D X = *S X (A.5) 

is the gauge invariant monopole creation operator. It creates the cloud of photons and the 
monopole at the point x. In ( |A.5| ) *5 X is the lattice 5-function, it equals to unity at the 
site x of the dual lattice and is zero at the other sites. Note that in the above formulas the 
radial part of the Higgs field which carries the magnetic charge is fixed to unity. 



Coming back to the original partition function ( A.l ) we get the expectation value of 



the monopole creation operator in terms of the fields 9: 



< $ m on{x) >=\\ V0exp{-S(d9 + W P )} , W P = 2tt5A-\D x - u x )), (A.6) 

— 7T 

where the Dirac string attached to the monopole ||, is represented by the integer valued 
1-form *uj x , which satisfies the equation: S*uj x = *5 X . 
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Appendix B 



Here we represent the partition function (A.l) for the compact electrodynamics with the 



general type of the action as the sum over the monopole currents. First we insert the 

+00 _ ^ 

unity 1 = J T>G 5(G — n) into the sum flOD and use the Poisson summation formula: 



E„ S(G - n) = En e 2wi ( G ' n \ We get: 

Z = const -Jv9 J VG F(G) eyip{i(d6 + 2irn,G)} . (B.l) 

Here G is a real-valued two-form. It is possible to change the summation variable n: 
/J/( n ) = + d<?), where n = m[j] + dg , dm [7] = j , dj = 0. Now we 

n q 5*j=0 

E+vr 
J T>9 f(d9 + 

n _7r 

+00 

27m) = ^ / T>A f(dA + 2n8A~ 1 j), where A = 9 + 27rA~ 1 5m[j] + 2irq and we use the 

<5*i=0-°° 

Hodge-de-Rahm formula: m = 5A _1 j + dA~ 1 5m. The integral over A gives the constraint 
5(5G), which we solve introducing the new integration variable H, G = 5H. Taking into 
account the relation d5A~ 1 j = j, valid for any j, such that dj = 0, we finally get the 
representation of the partition function as a sum over the conserved monopole currents: 

Z = const- Y e- Smon{ * j) . (B.2) 



*i(* Cl )e^ 



where 



Smon(*j) = - ln i j VHF{5H) exp {27ri(*H, . (B.3) 

The monopole action is nonlocal due to the integral over H. If we start from the 
Villain action S v {dO) = — In E n ex P { — @\\d9 + 27rn|| 2 }, then the integral over H in 
( |B.3| ) is Gaussian, and we get the well known expression for the monopole action: 
Smon{*j) = 47r 2 /3(*j, A~ l *j). Using transformations similar to (|B.l|)-(pl2D for the expec- 
tation value of the monopole creation operator (^), we get 

< $mon(x) >= — Y exp{5 m0 n(*j - * D x )} , (B.4) 

$*j=Sx 



where D x is defined by eq. ( |A.5| ). Therefore $ mon (^) creates a non-closed monopole 
world trajectory starting at the point x, which shows that, indeed, &mon(%) is a monopole 
creation operator. 
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Appendix C 



Let us briefly summarize the main notions from the theory of differential forms on the 
lattice [ p~9|| . The advantage of the calculus of differential forms consists in the general 
character of the expressions obtained. Most of the transformations depend neither on 
the space-time dimension, nor on the rank of the fields. With minor modifications, the 
transformations are valid for lattices of any form (triangular, hypercubic, random, etc). A 
differential form of rank k on the lattice is a function <^ defined on fc-dimensional cells 
of the lattice, e.g., the scalar (gauge) field is a 0-form (1-form). The exterior differential 
operator d is defined as follows: 



(d0)( Cfc+1 )= £ <K Cfc ). (c.i) 

Here <9c& is the oriented boundary of the fc-cell c^. Thus the operator d increases the rank of 
the form by unity; dip is the link variable constructed, as usual, in terms of the site angles <p, 
and dA is the plaquette variable constructed from the link variables A. The scalar product 
is defined in the standard way: if cp and ifi are fc-forms, then ((p,ifi) = J2c k ( / 3 ( c /c)V'( c fc)) 
where J2 Ck ^ s the sum over all cells c^. To any fc-form on the D-dimensional lattice there 
corresponds a {D — /c)-form *$(*Cfc) on the dual lattice, being the (D — fc)-dimensional 
cell on the dual lattice. The co-differential 5 = *d* satisfies the partial integration rule: 
(ip, Sip) = {dip, if)). Note that <5$(cfc) is a (k — l)-form and S^(cq) = 0. The norm is defined 
by: ||a|| 2 = (a, a); therefore, \\B — dip + 2irn\\ 2 in ( A.4 ) implies summation over all links. 
J2l( Cl )e% denotes the sum over all configurations of the integers I attached to the links c\. 
The action ( A.4 ) is invariant under the gauge transformations B' = B + da, ip' = ip + a due 
to the well known property d 2 = 5 2 = 0. The lattice Laplacian is defined by: A = d5 + 5d. 
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Figure 1: F($) for the confinement (a) and the deconfinement (b) phases. 
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Figure 2: The dependence of <& c on the spatial size of the lattice for three values of (3. 
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Figure 3: The dependence of 3>* n ^ on j3. 
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